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LONG-TERM TEMPERATURE EFFECTS ON GaAs SOLAR CELLS 


By 

John H. Heinbockel 1 and Ki H. Hong 2 


SUMMARY 

The thermal degradation of AlGaAs-GaAs solar- cells resulting from 
a long-term operation in a space environment is investigated. The solar 
cell degradation effects caused by zinc and aluminum diffusion as well 
as deterioration by arsenic evaporation are presented. 

Also in this report, we present the results of experimental testing 
and measurements of various GaAs solar cell properties while the solar 
cell was operating in the temperature range of 27° C to 350° C. In 
particular, the properties of light current voltage curves, dark current 
voltage curves, and spectral response characteristics are presented. 

Finally, s.ome theoretical models for the annealing of radiation 
damage over various times and temperatures are presented. 

INTRODUCTION AND STATEMENT OF PROBLEM 

It is the purpose of this report to obtain fundamental information 
on the long-term effects of high operating temperatures on GaAs solar 
cell characteristics. Research was conducted in order to identify and 
obtain limits for the magnitudes of probable high-temperature effects 
on the GaAs solar cell characteristics. 

In order to summarize the thermal degradation effects on GaAs solar 
cells: (1) a literature survey was performed, (2) basic laboratory experi- 

ments were performed to obtain the necessary fundamental information that 

1 Professor, Department of Mathematical and Computing Sciences, Old Dominion 
University, Norfolk, Virginia 23508. 

2 Research Professor, Department of Mathematical and Computing Sciences, Old 
Dominion University, Norfolk, Virginia 23508. 



will be needed for the design and usage of GaAs solar cells in a space 
environment over a 10-year time period, and (3) theoretical models were 
developed in order to estimate the removal of radiation damage effects 
due to annealing at high temperatures. 

Discussion is divided into three parts concerned with (1) thermal 
degradation of AlGaAs-GaAs solar cells due to zinc and aluminum diffusion, 
(2) experimental results of testing AlGaAs-GaAs solar cells at high 
temperatures (27° to 350° C) , and. (3) development of theoretical models 
for the prediction of annealing effects on the removal of radiation damage 
experienced by GaAs solar cells operating in a space environment. 
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THERMAL DEGRADATION OF AlGaAs-GaAs SOLAR CELLS 


Recently many reports on radiation damage of AlGaAs-GaAs solar cells 
in a simulated space environment have appeared (refs. 1 to 3). We are 
concerned about the thermal degradation of GaAs solar cells over a 10-year 
period of operation in a space environment. During this project, we studied 
the thermal deterioration of AlGaAs-GaAs solar cells which could result 
from a long-term operation in space. In order to assess the potential 
severity of the solar cell degradation, we assumed that the solar cell 
operated at a temperature of either 250° C or 350° C and in high vacuum 
conditions, i.e., 10" 10 to 10" 11 torr. 

The AlGaAs-GaAs solar cell has a three- layer structure. The bottom 
layer is nGaAs which is 200 to 350-pm thick; The middle layer contains 
zinc doped pGaAs and is approximately 1-pm thick. The top layer is a 
conducting layer of pAlGaAs grown by the liquid phase epitaxy (LPE) method 
and doped with zinc and is less than 1-pm thick. 

We investigated three potential mechanisms which may cause the solar 
cell to permanently deteriorate. These are (1) zinc diffusion for the 
Zn-doped p- layer into nGaAs, (2) aluminum diffusion from the LPE grown 
layer into the substrate, and (3) arsenic evaporation. For this study 
we employed currently available empirical data to predict potential deterior- 
ation. For the Zn diffusion we assumed that the total thickness of the Zn- 
doped p-layer is about 1 pm. Zinc diffusion during operation at temperatures 
in the range of 250° to 350° C for a period of 10 years was expected to 
increase the junction depth. There are many reports on Zn diffusion into 
GaAs (refs. 4 to 6) . In this study we selected two different Zn diffusion 
coefficients. These diffusion coefficients were obtained from curve fitting 
of data from both elemental (ref. 4) and ternary sources (ref. 5). One 
diffusion coefficient was for temperatures below 800° C and the other was 
for temperatures below 744° C. Both data curves were extrapolated to 
350° C and 250° C in order to estimate the junction depth change AX . , 
which was assumed to be the distance from a step function gradient in 
the zinc concentration. 'The resulting junction depth changes, which are 
given in table 1, agree within 30 percent of each other. We believe that 
this 30 percent difference results from the experimental error or the 



inaccurate determination of the diffusion coefficients during the curve 
fitting. The changes depicted in table 1, of approximately 0.1 pm, are 
small compared with the total p- layer thickness. However, this change of 
0.1 pm could be important for radiation resistance cells having junction 
depths less than 0.5 pro. 


Table 1. Thermal degradation of AlGaAs-GaAs solar cell for 
10 years of operation. 



Change in Junction Depth 

Arsenic Evaporation 

Temperature 

AXj 

(mg) 


C°c) 

(pm) 

GaAs 

Zn-GaAs 

250 

6.98(10' 4 ) a 8.81 (10 -4 ) b 



350 

5.78(10~ 2 ) a 9.9 (10" 2 ) b 

4.45(10~ 9 ) C 

6. 76(10 _7 ) C 



No arsenic left** 



Reference 4. ^Reference 5. Reference 8. ^Reference 9. 


For aluminum diffusion there is the problem of the movement of the 
GaAlAs-GaAs interface. We assume that the top layer of an AlGaAs-GaAs 
solar cell is composed of (Al^Ga^ with y = 0.9. Chang and Koma (ref. 7) 
have recently studied interdiffusion of A1 and Ga in (Aly.Gaj_y.)As as a 
function of temperature and aluminum composition. The aluminum diffusion 
coefficient and activation energy in GaAs can be compared with the zinc 
diffusion coefficient in GaAs. The thermal diffusion coefficient (D A1 ) 
of aluminum in GaAs at 350° C is approximately 1.7(10 -31 ) cm 2 /sec, whereas 
the thermal diffusion coefficient of zinc (D^) in GaAs is approximately 
10" 20 cm 2 /sec. The activation energies for diffusion of Al, Ga, and Zn 
in GaAs are 3.6 eV, 4.2 eV, and 2.49 eV, respectively. The aluminum 
diffusion coefficient is approximately 10“ 11 smaller than the zinc 
diffusion coefficient; therefore, a junction depth change due to aluminum 
diffusion can be ignored. 



Arsenic evaporation can occur when an AlGaAs-GaAs solar cell is 
operated at high temperatures. This is because arsenic vapor- evaporates 
from the surface of the solar cell at high temperatures, and this decompos 
the GaAs. The arsenic loss from a GaAs surface creates crystal defects, 
and these defects form recombination centers and alter the composition of 
the GaAs surface structure. The arsenic contained in a 1-cm 2 , 1-ym thick 
GaAs layer is approximately 0.7 mg. Based. on recent kinetic studies, we 
can calculate total arsenic loss for AlGaAs-GaAs solar cell when the solar 
cell operates at 350° C in space for 10 years. Using the results of Lou 
et al. [ref. 8), we estimate the total arsenic loss to be approximately 
6.7 x 10“ 7 mg over 10 years. Accordingly the density of arsenic vacancy 
is approximately 5 x l0 16 /cc. This density is a factor of 10 2 lower than 
the dopant density in the p-layer, and hence the total arsenic loss is 
negligible and a surface stoichiometry is maintained. On the other hand, 
total arsenic loss based on Arthur's results (ref. 9) suggests a more 
serious problem. The estimated arsenic loss rate from the GaAs surface 
is 1.25 x io-7 mg/cm 2 -sec, which implies that all arsenic in a l-pm layer 
of GaAs could evaporate within 2.5 months. Recently Sebestyen et al. (ref 
10) reported the arsenic loss from metal- covered GaAs. The calculated 
arsenic loss based on these experimental arsenic evaporation rates is 
shown in table 1. Based on the results of Lou et al. and Sebestyen et al. 
(refs . 8 to 1-0) , it is difficult to have confidence in any estimate of 
arsenic evaporation rates. 

As a direct observation of solar cell deterioration, we heated an 
AlGaAs-GaAs solar cell to approximately 300° C in a gas-tight enclosure 
with a quartz window for approximately 100 hours. From our direct obser- 
vation of photovoltaic parameters, there was no significant degradation 
in the solar cell characteristics due to short-term exposure to high 
temperature. Consequently, we feel that further work is needed on the 
problem of arsenic evaporation in order to quantitatively predict long- 
term arsenic loss (surface defects on GaAs and composition change) at 
high temperature. u 



THERMAL PROPERTIES OF AlGaAs-pGaAs-nGaAs SOLAR CELLS 


Introduction 

Gallium arsenide (GaAs) solar cells are being considered for use in 
generating electrical power in space on a solar power satellite. They 
are also being considered for use In a high- temperature, near-Sun space 
science mission. The space environment in which the solar cells will be 
operating will have high temperatures and high doses of radiation which, 
over long periods of time, degrade the operation of conventional solar 
cells. Gallium arsenide solar cells have certain advantages over conven- 
tional solar cells. These advantages are (1) a higher band gap, (2) higher 
open circuit voltage, (3) capability of operating at a higher temperature, 

(4) a high resistance to radiation, and (5) a spectral response which more 
closely matches the Sun's solar spectrum. 

In the sections that follow we will be concerned with the operational 
characteristics of pAlGaAs-pGaAs-nGaAs solar cells operating within the 
temperature range of 27° C to 350° C. Test data was obtained from GaAs 
solar cells operating in a simulated space environment. This test data 
characterizes the photovoltaic parameters which, as functions of temperature 
determine the behavior of GaAs solar cells. 

The test data, recorded as functions 'Of temperature, consisted of 
illuminated current voltage curves, dark current voltage curves, and 
spectral response curves. 

Experimental Considerations 

The following describes the experimental work that was performed in 
order to determine the high-temperature GaAs solar cell thermal properties. 
Ten heteroface pGaAlAs/pGaAs/nGaAs solar cells were fabricated at NASA/ 
Langley Research Center using the etch-back epitaxy process. The n-type 
GaAs substrates were Si doped to a carrier concentration of 2 x 10 17 
carriers/cm 3 . The GaAlAs layer thickness was between 0.3 and 0.5 pm. 

The pGaAs layer was approximately 1-pm thick, giving a junction depth 
of 1 pm (as measured from the GaAlAs interface) . Large area, vacuum- 
evaporated, Sn-Ag contacts were used for electrical connections to the 
nGaAs. Front finger contacts were deposited to the GaAlAs layer by 



sputtering a layer of Pa followed by a layer of Ag using an appropriate 
mask. 

Before operating the solar cells at a high temperature, two GaAlAs- 
GaAs solar cells were prepared in order to test the stability of the solar 
cells upon heating. The purpose of the stability check was to obtain the 
open circuit voltage change due to the current leakage during the heat 
treatment procedure. There was no observable change in the open circuit 
voltage, short circuit current, or fill factor. 

The solar cells were mounted on stainless steel headers and were 
given an antireflection coating. Before they were heat treated, their 
photovoltaic properties were measured and recorded. Heat treatment 
studies at 200° C were conducted by separating the solar cells from the 
headers and placing them in a vacuum oven for approximately 5 hours. 

During the heat treatment of these solar cells, the temperature was 
continuously recorded. After heat treatment the samples were remounted 
on the stainless steel headers and the photovoltaic properties were 
recorded. In particular, the open circuit voltage, short circuit current, 
fill factor, efficiency, spectral responses, and dark current voltage were 
measured both before and after heat treatment. The cells tested proved 
to be stable. 

For the study of high- temperature operation of GaAs solar cells, 
a small electric heater together with a stainless steel header mounting 
plate and chromel-alumel thermocouple wires were prepared and installed in 
a. small, blackened;-.' stainless steel cylindrical chamber. Desired tempera- 
tures were obtained by adjusting the current input with two variacs. 

Alumina and titanium oxide were used to blacken the stainless steel 
cylindrical chamber in order to reduce reflection of light. The proper 
electrical insulation was achieved by using AI 2 O 3 . During the measurements, 
the solar cells were placed in a stainless cylindrical chamber with a 
quartz glass window. which allowed the solar cells to be irradiated. A 
thermocouple wire was attached underneath to the stainless steel header 
plate in order to monitor the temperature during each measurement. Dry 
argon was slowly flushed into the system in order to maintain a small 
positive pressure to prevent the solar cell junction from oxidation. 

The temperature measurements on the solar cell were carefully calibrated 



with special equipment designed in the laboratory. The temperature 
stability was maintained with less than ±1° C during the measurements. 

Since silver epoxy was not adequate for the solar cell operation 
at the temperatures above 180° C, the thermal compression bonding tech- 
nique was employed for sample mounting and electrical contacts. Usually 
3 to 4 gold wires having a 2-mil diameter were used to attach the cell to 
the header. Before mounting the cell on the header, the headers were 
coated with a thin gold film for better adhesion. The heat conduction 
through the gold wire to the solar cell was carefully observed by moni- 
toring the temperature on the header and on the solar cell junction. 

During all measurements, sufficient time was allowed for the solar cell 
to establish itself at an equilibrium temperature. 

Measurements were repeated and the reproducibility of the measure- 
ments was maintained. For the -illuminated current .density voltage measure- 
ments, a solar simulator light source was used to produce an input power 
of 135.3 mW/cm 2 at AMO. For the dark current density voltage measurements, 

the V -J method was employed in order to reduce the effect of the series 

oc SC c j 

resistance. 

The spectral responses as a function of temperature were measured 
using a mercury lamp, electrometer, and a set of narrow band pass filters. 
The mercury lamp was chosen in order to obtain a strong light intensity 
for all wavelengths. This also helped to compensate for the increasing 
dark current with increasing temperature and thereby produce enough 
current density at high temperatures for measurement. To achieve mono- 
chromatic light, a series of narrow band pass filters with approximately 

O 

50 A resolution was used. The current and photon density were recorded 
at each wavelength, enabling us to obtain absolute spectral responses. 
Measurements were repeated several times for each solar cell in order 
to ensure maximum accuracy. The wavelengths chosen for the spectral 
responses were from 0.44 ym to 1 ym with 0.04 to 0.05-ym v/avelength 
intervals and temperatures ranging from 28° G to 346° C with approxi- 
mately 60° C incremental steps. The absolute spectral responses were 
numerically integrated to produce the current density J sc , which was 
then compared with the J obtained from the illuminated J-V measure- 
ments. The obtained values for the current densities were well within 


9 



the error range involved in the approximation. All the measurements were 
repeated, and reproducibility was carefully checked at each time of the 
run. 


Light Current Voltage Curves 

The current voltage (J-V) characteristics under AMO sunlight for 
GaAlAs-GaAs solar cells in the temperature range from 27° to 350° C were 
studied. Similar measurements were made by Hovel (ref. 11) by employing 
a xenon light source in the temperature range from 27° to 245° C. In 
the past, the values of illuminated photovoltaic characteristics above 
T = 245° C were obtained by extrapolation. This was because of the 
experimental difficulties encountered at temperatures above T = 245° C 
(refs. 11, 12). In this study we extended the temperatures beyond T = 
245° C and made actual measurements of illuminated photovoltaic charac- 
teristics of GaAs solar cells. up to a temperature of T = 350° C. For 
these measurements the solar simulator was used to reduce the excessive 
intensity of light at 0.82 ym. The measurements obtained were not 
corrected for the contact area masking. 

Figure 1 summarizes the light current voltage curves for temperature: 
in the range from 27° to 350° C.- The rest of the cells tested exhibited 
a similar behavior. The following are some nominal values recorded at 

27° C: V qc = 0.97 V, FF = 0-.7, J sc = 24 mA/cm 2 , and efficiency = 

n = 12 to 14 percent. Using TiC >2 antireflection coating, the J gc 
increased to 33 mA/cm 2 and the efficiency n rose to nearly 16 percent. 
As the temperature of the solar cell increased, the open circuit voltage 
decreased due to the dark current increase. The rate of change for V qc 
over the entire temperature range was observed to be AV/AT = -2.2 mV/°C 
for all solar cells tested. This value is in agreement with previous 
measurements (ref. 11). We also observed that the short circuit current 
density increased with increasing temperature, but not in the manner 
reported earlier (ref. 11) . The short circuit current density changed 
with increasing temperature in the following manner: the short circuit 

current 'density increased with increasing temperature due to the improve- 
ment in the diffusion length and a shift in the absorption edge to lower 
energies. The short circuit currents at 27° C for all solar cells were 



Figure 1. Photovoltaic parameters as a function of temperature. 



observed in the neighborhood of 24 mA/ cm 2 . The short circuit current 
gradually increased with increasing temperature up to around 210° C. Then, 
for temperatures above 210° C, the short circuit current density slowly 
decreased with increasing temperature. At temperatures above 270° C, there 
was a considerable decrease in the short circuit density. Figure 2 demon- 
strates that the spectral response at 346° C was reduced. This reduction 
could be attributed to an increased dark current or to contact resistance. 

In general our measurements for the change in short circuit current density 
with temperature AJ /AT showed much smaller rates than that obtained in 

SC 

reference 11. In the temperature range 27° C to 210°' C we observed that 
AJ sc /AT = 2 mA/cm 2 . 

For most solar cells the fill factor at 27° C was-- in- the neighborhood 
of 0.65 to 0.7, as compared to 0.78 in reference 11. As temperature increased, 
the fill factor decreased, and at T = 350° C we found the fill factor to be 
around 0.35. 

The efficiency at 27° C was between 12 and .14- percent for most of the 
solar cells tested. We observed that the efficiency changed with tempera- 
ture and decreased at the rate of 0.035 percent/°C. This value is slightly 
higher than the rates of 0.02 to 0.03 percent/ fl C reported in reference 11. 

For our measurements, the small improvement in the short circuit current 
density did not offset the linear decrease of the open circuit voltage 
with increasing temperature. Accordingly, the efficiency decreased 
linearly with increasing temperature. The measured efficiency at 350° G 
for most solar cells with antireflection coatings was in the neighborhood 
of 2 percent. The GaAlAs-GaAs solar cells were usable at and up to the 
temperature of T = 350° C, but the measured values of efficiency were 
lower than the predicted values obtained from the extrapolation of - 
earlier measurements (refs. 11, 12]. 

Dark Current Voltage Curves 

The dark currents in GaAs solar cells are composed of injected 
majority currents over the potential barrier and the recombination 
current of holes and electrons within the depletion region. The dark 
current voltage characteristics are determined by the combined effects , 
of the current transport mechanisms and series and shunt resistance. 
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SPECTRAL RESPONSE 
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Figure 2. Measured spectral responses as a function of temperature. 


The details of the dark J-V mechanisms are well described elsewhere . 

(ref. 13). The injection currents vary as exp[qV/KT], and the recombina- 
tion currents vary as exp [qV/2KT] . Accordingly, the recombination currents 
dominate at a low forward bias, and the injection current dominates at a 
higher bias. For a high band gap material such as GaAs, the recombination 
current is generally greater than the injection current. As the tempera- 
ture increases, the intrinsic carrier density increases due to the band gap 
decrease. The dark current largely determines the voltage output. Hence, 
with increasing temperatures, the output voltage decreases due to the 
increasing intrinsic carrier densities. 

s. 

The data collected was obtained using the V -I method. This- method 

6 oc sc 

minimizes the effect of the series resistance on the measured currents and 
voltages. The current voltage data as a function of temperature is illus- 
trated in figure 3. Most measurements show two exponential regions exist 
which indicate the two different current components. The slope of the J-V 
curves determines the diode A-factor and saturation current density J q 
at zero bias. An A-factor of A = 2.1 was obtained for most solar cells at 
room temperature. In general, the A-factor obtained from the experimental 
data was slightly higher than the predicted value. Some solar cells had 
an A-factor as high as A = 2.6. This high an A-factor indicates the 
existence of poor diode characteristics. Metal precipitation, oxygen, 
and native defects are partially responsible for this high an A-factor 
value. During this study, the A-factor and J q values were determined 
by single exponential least squares curve fitting of the data, for which 
the mean square difference between the measured and computed values was 
minimized. The A-factor values and J Q values thus obtained are given 
in figure 3. The A-factor values decrease from 2 to 1.6 with increasing 
temperature. This is to be expected since the ideal diode A-factor is 
expressed as 


A 



1 

AlnJ 

AV 


where A is the difference between two data points. 


CD 
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Since the data indicates that the recombination current is dominant, 
a composite lifetime t c as a function of temperature can be obtained 
(ref. 14) as 


where J q is the saturation current density, tu is the intrinsic carrier 
density, and W is the width of the depletion region. The depletion 
region can be obtained from the known doping density. The temperature- 
dependent value of W was calculated assuming a slowly varying dielectric 
constant. The depletion region is approximately 0.1 ym at T = 27° C. 

With increasing temperature the depletion region decreases to 0.08 ym. 

Using equation (2) , the composite lifetime at 18° C is approximately 
2(10 -9 )sec. Summarized in table 2 are the composite lifetimes (as a 
function of temperature) deduced from the dark J-V curves. 


Table 2. Recombination lifetime deduced from 

dark J-V as a function of temperature. 


Temperature 

Lifetime 

C°C) 

'(sec) 

18 

2.31 

X 

10~ 9 

90= 

1.75 

X 

10-8 

150 

5.9 

X 

10~ 8 

210 

7.4 

X 

10' 8 

270 

1.04 

X 

10“ 7 

350 

1.55 

X 

O 

1 


• This composite lifetime (ref. 14) is a weighted average of the minority 
carrier lifetimes in the p and n materials. In the next section these 
composite lifetimes will be compared with lifetimes obtained by fitting the 
spectral response data. 



• Spectral Response Curves 

The spectral responses o£ the GaAs solar cells were measured at 
different wavelengths and different temperatures. The experimental data 
obtained was read by an HP9830 calculator to produce the measured absolute 
spectral responses. The measured spectral responses were observed to be 
flat between the 0.55 ym and 0.85 pm wavelengths for all temperatures 
considered. 

In order to ensure that we had absolute spectral response measurements, 
we numerically integrated the spectral response curve with the solar spectrum 
over the entire wavelength in order to obtain the current density J 

S 0 

The value of J sc obtained in this way agreed within 10 percent with the 
current density obtained from the illuminated current voltage data. The 
small discrepancy between the measurements was probably caused by the 
approximate integration procedure used over the entire wavelength. 

The change in the spectral response at the blue end of the spectrum, 
with increasing temperature, was small when compared to the spectral 
response change at the red end of the spectrum. The spectral response 
in the intermediate wavelength regions lying between the blue and red 
regions was essentially a constant. The current, in the red end of the 
spectrum, increased with increasing temperature due to the band gap 
reduction with increasing temperature. However, the spectral response 
with increasing temperature did not significantly affect the total 
current. These measurements give the complete spectral response as a 
function of temperature for the first time. 

In order to explain the measured spectral response curves with 
increasing temperature, we used the Hovel model for the spectral response 
as developed in reference 15. In this analytical model for the spectral 
response we made many assumptions. Since the absorption coefficient of 
AlGaAs as a function of temperature is not known, we assumed that this 
absorption coefficient (along with the absorption coefficient of GaAs) 
behaved like the band gap of GaAs as a function of temperature. The 
following temperature relationships - were also employed- in the model. 
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D a T -1 / 2 , u a T _3/ ^ 2 

P 

T a T 3 / 2 , a T -3 / 2 C3) 

L a T^/2 
/ 

where D is the diffusion coefficient, t is the lifetime, L is the 
diffusion length, T 'is the temperature in degrees Kelvin, and p p , u n 
are mobilities. 

The calculated spectral responses were compared with the measured 
spectral responses. The best fits of the experimental spectral responses 
as a function of temperature are illustrated in figures 3 and 4. The 
parameters used -in the best fit model are summarized in table 3. The 
parameters that are critical in fitting the blue end of the spectrum are 
the aluminum composition, the recombination velocity, the diffusion length, 
and lifetime in the AlGaAs . At the red end of the spectrum, the parameters 
most important are the hole diffusion length and lifetime. In general, 
there was good agreement between measured and calculated spectral responses. 

The band gap of AlGaAs as a function of temperature as well as aluminum 
content in AlGaAs and surface recombination velocities in AlGaAs and pGaAs 
must be accurately determined to further reduce the error in fitting the 
spectral response curves at the blue and red ends of the spectrum. 

Various solar cell parameters were estimated in obtaining the computer 
fits to the experimental data. These values are given in table 3. Note 
that the recombination lifetime was also obtained in this way from the 
spectral response measurements. The value obtained for the lifetime by 
this method differed by a factor of 10 from t-he composite lifetimes 
calculated earlier. 

THEORETICAL MODELS FOR HIGH-TEMPERATURE OPERATION OF 

Ga, A1 As-GaAs SOLAR CELLS 
1-x x 

In this section we present some theoretical models for the high- 
temperature operation of GaAs solar cells. This section is based upon 
the discussions found in reference 16 through 20. 
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Figure 4. Spectral response curves at various temperatures 
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Table 3, Parameters deduced from temperature-modeled spectral responses. 


Temperature t (sec) D t (sec) D t (sec) D 

(°K) ' L a (cm) (cm 2 /volt-sec) L g (cm) (cm 2 /volt-sec) l p (cm) (cm 2 /volt-sec) 

a ■ g p 



1 

X 

10- 7 



IO " 3 

1.13 

X 

o 

I 

H -* 

o 

88.70 

4. 66 

X 

10 - 9 

7.75 

300 



6.25 

X 




2.5 

X 

10 - 5 




1.0 

X 

10- 4 


1.90 

X 

10 -4 


363 

1 

:< 

10- 7 

6.25 


o 

1.50 

X 

10-10 

80:62 

6.20 

X 

10-9 

7.04 



X 

IO - 3 




2.5 

X 

10" 5 




1.1 

X 

io - 4 


2.09 

X 

10- 4 


423 

1 

X 

10" 7 

6.25 



1.89 

X 

o 

r — 1 
1 

o 

I — 1 

74.68 

7.8 

X 

10“ 9 

6.52 


* 

X 

10 3 





2.5 

X 

io - 5 




1.19 

X 

10“ 4 


2.26 

X 

10 " 4 


481 

1 

X 

io - 7 

6.25 

X 

10 “ 3 

2.29 

X 

10" 10 

70.04 

9.95 

X 

io -9 

6.12 


2.5 

X 

10 " 5 




1.27 

X 

10~ 4 


2.40 

X 

10" 4 


540 

1 

X 

10" 7 

6.25 

X 

f — ' 

O 

l 

CO 

2.72 

X 

10 - 1 ° 

66.099 

1.13 

X 

10-8 

5.78 


2.5 

X 

10- 5 




1.34 

X 

IO" 4 


2.55 

X 

IO" 4 


619 

1 

X 

IO -7 

6.25 

X 

10-3 

3.34 

X 

10- 10 

61.74 

1.38 

X 

10“ 8 • 

5.39 


2.5 

X 

10- 5 




1.44 

X 

10- 4 


2.73 

X 

10" 4 



Fixed Parameters Used for Calculation 

S & = 10 6 cm/sec = 10 4 cm/sec 

Hole density = 2 x 10 8 . Electron density ~ 1.3 x 10 17 

Window thickness = 0.3 pm Junction depth = 1 pm 

x in (Al Y Ga 1 _ v )As = 0.95 


to 

cn 



The interaction of electron irradiation with the GaAs solar cell pro- 
duced displacements of atoms within the crystal lattice. The fraction of 
these displacements which remain after long periods of time at room tempera- 
ture is called "permanent damage." The process of annealing point defects 
involves their diffusion in the lattice until they encounter traps or sinks 
where they will then cease their meandering within the crystal. 

In the electron irradiation of semiconductor crystals, coulomb scatter- 
ing events will occur and the electrons will transfer some of their energy 
to the lattice atoms. The maximum amount of energy that can be transferred 
is given by (ref. 16) 

2[E + 2m c 1 2 )E 

T = e e e 

111 Me 2 

where E = electron energy, m = electron mass, c = velocity of light, 

C 6 

M = lattice atom mass. For a GaAs crystal the maximum energy that can 
be transferred by 1-MeV electrons is 


T _ 2[1 * 2(0.510577)] CD, 
V 64890.442 


62.294 (eV) 


T ,As = 2 
m* 


1 + 2(0.510577) 

69731.724 


(1) = 57.969 (eV) 


Due to long-range interaction of coulomb collisions, the average amount of 
energy transferred to the lattice is less than these theoretical maximum 
values. There is a certain displacement threshold energy E^ which will 
displace a lattice atom to an interstitial position. This energy is between 
10 and 60 eV for most semiconductor materials. The displaced atoms will 
leave behind a vacancy in the crystal lattice, and the resulting inter- 
stitial vacancy pair is called a "Frenkel defect." For most semiconductors 
these defects are produced uniformly throughout the crystal lattice. The 
total number of defects produced by an electron irradiation is a function 
of the electron energy and can be expressed as 
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N = o d U^ 


( 6 ) 


where 

N = the total number of defects/cm 3 

= the effective displacement cross section (cm 2 ), which is 
a function of electron energy (see fig. 5) 

= the number of lattice atoms/cm 3 
<j> = the electron fluence (e/cm 2 ) 

As defects are produced, various electrical characteristics and 
responses of the GaAs solar cell are degraded. Annealing is the process 
of raising the temperature of the solar cell in order to impart an 
activation energy to the defects so that they become mobile within the 
crystal and will meander through the lattice seeking to obtain a lower 
energy state. Defects are removed by this annealing, and electrical 
properties that were degraded by the creation of defects are restored by 
the removal of defects. The amount of restoration in the electrical prop- 
erties depends upon the annealing time and the annealing temperature. 

The annealing of GaAs is essentially in two stages: (1) The Ga atoms 
require a fairly low activation energy which is around room temperature 
(27° to 127° C) , and at this annealing temperature the int-erstitial- 
vaaancy pairs combine and destroy each other. (2) The remaining damage 
is to the arsenic atoms which are displaced. These atoms require annealing 
around 240° C (513° K) in order to receive the activation energy necessary 
to become mobile. The annealing process at temperatures greater than 
240° C is fairly complex as can be illustrated by considering the various 
reactions that can occur during the different stages of annealing. 

Tlfr-x 1^4- 


Vj = vacancy 

V 2 = di-vacancy 

V = cluster of n-vacancies 
n 

il = interstitial 
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i = cluster of n- interstitials 
n 


Kj = reaction constants (rates) , 


then the annealing kinetics can be categorized by various types of reactions. 
Some of these reactions are 


(a) Recombination reactions. - 


'Ki 

(1) Vj + ij -> annihilation 

(recombination of closely spaced vacancies and interstitials) 


K 2 

(2) ij ->• sinks 

K 2 

(3) ii + ii 5 i 2 

t 3 

(4) ii + impurity j complex 

£5 

(interstitials combining with impurities) 


(5) higher order reactions. 


(b) Reactions due to diffusion of defects within lattice. - 


k 6 

(1) V X + Vi ± v 2 
Ki 


(2) Vi sinks 


(3) V 2 3- sinks 
k 8 


k io 

(4) Vi + impurity complex 

K tl 

(impurity atom can capture vacancies) 


(5) higher order reactions . 
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( c ) Clustering reactions. - In addition to the reactions stated in (a) 

" arid (b) , the presence of impurities will suppress various annealing 

reactions. These higher order and more complex types of reactions 
are' assumed to be secondary to the basic reaction of vacancies and 
interstitials annihilating each other. 

With reference to figure 6, let equal the number of defects produced 

by a 1-MeV electron fluence of <p (e/cm 2 ). Using equation (6) we have 


N D = 


w 


(7 


0 


► DEFECTS PRODUCED 

H 1 

N N 

r 

U DEFECTS REMOVED 



Figure 6. Defect production and removal. 


Using t as time, N = N(t) the defects at time t, T = temperature (°K), 
and as the defects which remain after a long time, we assume that the 

recovery rate is proportional to the defects remaining. This assumption is 
analogous to assuming a first order or monomolecular chemical reaction 
between the vacancies and the interstitials. We write 


jL(M D - N) a Mp e KT CM - up 

recovery proportionality defects 
rate constant remaining 


( 8 ) 


where 



N = number of defects at any time t of the annealing process 
N d = = defects produced by fluence <j> 

-E /KT 

3 . 

&> o e = Boltzmann factor with E & = activation energy (eV ) , 

K ■ is. Boltzmann! s -constant, T is annealing temperature 
(°K) , and cj^ is the frequency factor, which represents 
the frequency at which vacancies and interstitials 
annihilate each other 

N r = defects remaining after a long anneal time. 

The differential equation (8) is subject to the conditions N(0) = 

and limN(t) = N . It is readily verified that the solution of equation (8) 
t-*» 

which satisfies the above conditions is 

/ -E /KT\ 

N = N(t) = N r + (N d - N.^ exp|-<o o t e a -J (9) 

Let f denote the fraction of initial defects which is removed during 
an anneal time t. For = initial defects and N = N(x) = defects remaining 
after an anneal time x, we have 


f = 


N r 


N(t) 


N, 


CIO) 


or 


N = N(t) = N d (1 - f) 


Cll) 


How defects are removed depends upon the type of reactions that take 
place. . For a monomolecular reaction we have shown N = N(r) is given by 
equation (9) . For a bimolecular reaction, we can write 

|jr = F - aN - bN 2 (12) 

where F, a, b are constants and x is the annealing time. 
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Let 


(u = /a 2 + 4bf = a) exp{_E /KT> 
o a 

then we may express equation (12) in the form 



bdx 


( 13 ) 


An integration produces the solution 


N = N(t) 



(14) 


where N(0) = N D and A = INp + - u 


,N D + 


a + to 
2b . 


is a constant. This 


solution enables us to write equation (10) in the form 


f = 




(15) 


which represents the fraction of defects remaining in the case of biraolecular 
reactions. In the case where b tends to zero, then to = a and A = 0 so that 
for monomolecular reactions equation- (15) takes on the form 



where f = (Np - N r )/Np is the fraction of defects remaining after a long 
time. 

We now define as the equivalent fluence to produce N defects. 

Using equation (11) we may mite 


_ N(t) _ 1 


G d N * 


a d N * 


n d U 


f) 


(17) 
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where N = N(x) is a function of the annealing time x and represents the 
number of defects remaining after an anneal time x. 

The following is a simplified annealing model. We start with the 
approximate relation (ref. 15) 


I 

sc 


■*•300 


where F(x, S, 


F (x , S , L , X.) 
go go* go* y 

1 g g g r 

L, X) = sinh^-jr^ + cosh 



and 


( 18 ) 


go' 

go* 

go’ 


x are lifetimes (sec) before and after irradiation 
g 

S ct are recombination velocities (cm/sec) before and after 
irradiation 

L are diffusion lengths before and after irradiation (cm) 
g 

Xj is junction depth (cm) . 


We assume further that the lifetime, recombination velocity, and 
diffusion length change with fluence according to the relations 

1 1' 

— = — + D K $ 

x t g g e 

g go 5 s 


1 1 

— = + K $ 

L 1 2 * * L 2 * S 




go 


S = S , T 

g go\ L 


L N 2 
go 


g, 


(19) 


where D is the diffusion coefficient and K is a damage coefficient 
g g 

(ref. 15) . The quantity $ ^ is the effective fluence defined by equation 

(17) . Using equations (10) , (15) , and (17) we may write 


$ 


e 



( 20 ) 
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When equations (19) and (20) are substituted into equation (18), figures 

7 through 13 result. The nominal values employed to produce these graphs 

-E /KT 

are given. ;in figure 9. In these curves we have used w = io q exp( a' ) 

with (>) = 2.8766(10 12 ), f = 0.94. 

o v * ’ o 

As another way of predicting effects of annealing we will employ the 
following statistical models which were developed in reference 15 for 1-MeV 
electron fluences 4 and junction depths X ^ (pm) : 

075 logiof-M (21) 

VlO 16 / 


0.609 + 0.06 log 10 (— \ 

VLO 16 / 

( 22 ) 


logiol-r 2 *- 11 

sc 


= 0.8649 + 0.194Xj + 1, 


and 


'V 

•* i OCO q i 

10 SlO \-y U 

OC 


•0.6037 + 0.098X. + 
1 


x !og 10 




The current voltage relation for the solar cell is 

QV 


I = I - I 

SC C 


A KT 
o 


- 1 


(23) 


where when I = 0, V = V 


oc 


A KT 
o 


oc 


In 


sc 


= open circuit voltage and equation (23) gives 


+ 1 


(24) 


where A q is the A-factor (perfection factor), K is Boltzmann's constant, 
q is electron charge, I is short circuit current, and I Q is the dark 
current. To incorporate annealing effects, we replace § in equations (21) and 
(22) with the equivalent fluence defined by equation (20) . In addition. 
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FLUENCE (1 MeV ELECTRONS/ cm 2 ) 

Figure 7. Short circuit current vs. fluence before and after annealing at 200°C for r hours 
(theoretical curves A = 0) . 
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Figure 8. Short circuit current ratio vs. anneal time at 200°C (theoretical curves A = 0) . 
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Figure 11. Annealing of shallow junction (0.5-ym) cells at 200°C (theoretical curves A = 0) . 





we will let to Q = 0.077 exp{1.3/(8.6139(10“ 5 )T)} where T is the temperature 
in °K. From equation (24) we write 


sc 


I = 

O X 

e 0 - 1 


where x q - qV oc /A Q KT. Then equation (23) can be expressed as 


I = I - I 


x 1 
e - 1 


sc sc x 

i o i 

Ve - 1, 


where X = qv/A Q KT. This equation simplifies to 


I = 


sinhl 1 ? 


I Jl 

sc : 


-(X q -X)/2 


T77x \ 

sinh/ _ o 


(25) 


which is plotted in figure 14 for various values of V Q< J A ■ 
The cell power output is given by 


P = VI = VI 


sc 


rX 


1 - sinh 


V 


2 V 


OC/ 


sinhl 


\ 2 / 


("X 0 /2) Cl - V/V qc ) 


(26) 


dP 


where at maximum power P * P , I = I , V = V , and -jtj- = 0. Differen- 

max max in ax qv 

tiating equation (26) we obtain that at maximum power 


VI- r 
oc max _ g 

I V sinh ' g 

sc max 


exp 


V 


-di - 2- 


_max 1 

V 
, oc j 


(27) 


where 
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5 


X qV 

o _ oc 

2 2A KT 

o 


We may now use equation (25) to obtain the current voltage plot for irradia- 
tion for a particular cell/ and we can then use equations . (21) and (22) 
together with equation (25) to obtain after irradiation plots. These 
curves are illustrated in figures 15 through 18. 

The theoretical curves in figures 7 through 18 agree quite well with 
the experimental data obtained from various annealing experiments. 


Annealing During Irradiation 

Consider an in situ annealing experiment where samples are irradiated 
with 1-MeV electrons at a rate <ji (e/cm 2 -hr) while being held at a tempera- 
ture T (°K) . We asstune that production and removal of defects are a 

a 

simple monomolecular process and write 


rate of change J _ J defects \ I defects removed 
of defects J - | produced/ - (because of annealing 


dN - I a „ 

5t ■ W - % nsr N 

\ a; 


(28) 


which is subject to the initial condition N(0) = 0. 

In the case | is a constant, equation (28) has the solution 


N = 


cr,N.$ T 1 

d l , -cox 

1-6 

03 L 


(29) 


~ E a /KT a 

where t is the anneal time and co = w e .We have that 

o 


limit N = 


w 


oo 


0) 
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Figure 15. Current voltage characteristics of 10* 6 1-MeV electron/ cm 2 
irradiated GaA/As solar cell after t hours of annealing 
at 200°C obtained from statistic model. 
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Figure 16. Short circuit current for deep junction (4-yra) cells vs. 
o' • anneal time (statistical model) . 
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which says that after a long period of in situ irradiation at a constant 
flux rate $, there will be a constant damage to the solar cell. If t 
is the annealing time and also the irradiation time during which the solar 
cell receives a fluence of <{> (e/cm^), then we can approximate $ by 

<f>/t and write 


limit N 
t-*» 




(lit 


This implies that for a constant fluence $, the limiting damage will be 
smaller the longer the irradiation time t, and theoretically this limiting 
damage will be small after long annealing times and low dose rates. 


In-Flight Experiment 

During an in-flight space experiment it is assumed that the diffusion 
length L changes according to the relation given in equation (19), where 
4> e is an equivalent fluence defined by 


$ 


e 


N 


°d N * 


(30) 


It is assumed that the defects produced are dependent upon the operating 

temperature T q , the type of radiation, and the energy of the radiation. 

With reference to figure 19, it is assumed that the solar cells are operating 

at a temperature T q over a time period 0 < t a. During this time period, 

the defects produced by the irradiation increase. During the time interval 

a < t < a + t, annealing at a temperature T occurs which will reduce the 

number of defects. It is assumed that this process will be a periodic 

process, and we desire to study the defects produced as a function of 

o, t, T , and T . 

’ * o’ a 

Let 


T = T (t) 


T , 0 £ t <_ a 

T,o<t<o + x 
a 5 — — 


(31) 
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Figure 19. Radiation damage and 


2g+t 2ct+2t 



2g+t 2<j+2t 


as a function of time. 



with T(t + cx + x) = T(t). We assume that 


dN _ j amount \ f amount \ 
dt ~ ^produced J \removed J 


dN 
dt = 


4> d N £ 4 - w o Nf(t) 


where 


* = flux (e/cm 2 -hr) (assumed constant for orbital mission) 
N = number of defects/cm 3 
N = number of lattice sites/cm 3 
= frequency factor 
= displacement cross section 


f(t) = exp | | “ time-varying Boltzmann term where E & is 

activation energy and T(t) is given by equation (31). 
The solution of equation (32) satisfycing N(0) = N q is given by 


N = N (t) 


1 

g(t) 





f 


g(s)ds 


( 32 ) 


(33) 


where 


g(t) = exp- 


let 


G (t) = exp 


ir 

IX 


f(u)du> , 0 <_ t _< (cr + x) = P 


then, .since f(t) is periodic with period P, we may write 


(34) 


(35) 


51 



I 


g(t + nP) 


= exp|u o 


no e 


-a 


+ to nx e 


-6 


G(t) 


where 


(36) 


E 

a 



o 


and 



Then 


(37) 


N(t + nP) = 


g(t + nP) 


" _ /»t+nP -i 

N o + a d N A'J g( s ) ds J 


0 < t < P 


(38) 


In the limit as n increases without bound, the first term in equation (38) 
approaches zero, and the second term in equation (35) gives 


limit 


r Tt+nP 


X 


g(s)ds 


g(t + nP) 




ff dV 


-E a /KT(t) 


= N(t + nP) 


(u e 
o 


(39) 


Thus, after a long period of operation, the equivalent fluence is 


, _ $ 

^e ” » exp{-E & /KT(t)} 

which varies between the minimum value of 


6 . . 

■ e,mrn 


u e 
o 


and maximum value of 


<!> 


_ $ 

e,max -B 

5 oj e 

o 
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where a) Q is the frequency factor. The resulting damage to the diffusion 

length varies between L . , determined by 
& nun’ 1 


— — = — + K § 

1^2 -t 2 o e,min 
min o 


and L , determined by 
max 3 


= -i- + K $ 

L 2 . L 2 g e,raax 
max o 


If the junction depth X. is such that X.' < L . < L , then there will be 

) 1 mm max* 

no observable damage to the electrical characteristics, as any minority 
carriers produced by ' photon: impingement will be guaranteed arriving at the 
junction. 


CONCLUSIONS 

In this report we studied heteroface GaAlAs-GaAs solar cell operation 
in the temperature range 25° to 350° C. The striking difference between 
our work and previous reports is that we raised the operating -temperature 
beyond 250° C. We looked into illuminated light J-V, dark-J-V curves, 
and comprehensive spectral response characteristics as a function of 
temperature. From the illuminated light current voltage curves, the slopes 
AV/AT, AJ/AT, AFF/AT, and Ap/AT were measured and compared with previous 
reports . 

The actual measurement of efficiency near 350° C was observed to be 
approximately 2 percent at AMO, which is somewhat lower than the extrapo- 
lated values obtained from other sources. From the dark current voltage 
measurements, we observed that the A- factor value decreases from A = 2.1 
to A = 1.5. with increasing temperature (up to 350° C) . Assuming that A = 2 
for all temperatures, we could calculate the composite recombination life- 
time as a function of temperature. The difficulty in defining the slope 
of the ln(J) - V curve forced us to assume that A = 2 in order to obtain 
the recombination lifetime. Consequently, we could only obtain a crude 
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approximation for the lifetime using this technique. The recombination 
lifetime was also obtained from the spectral response measurements. The 
resulting values of the recombination lifetimes obtained by these two 
methods gave results which differed by a factor of 10. It is believed 
that the discrepancy between the values calculated by these two methods 
was due to the crude approximation in defining the A-value slope in the 
dark current voltage curves. 

We also obtained a comprehensive spectral response characteristic. as 
a function of temperature. The values of the spectral response were obtained 
both theoretically and. experimentally. In order to calculate the theoretical 
spectral response curves, the absorption coefficients of AlGaAs were modeled 
to change with temperature in a manner which agreed with the band gap change 
with temperature. A computer fit of experimental and calculated spectral 
response data produced many photovoltaic parameters as a function of .tempera- 
ture. We found that the aluminum content, the recombination velocity at the 
AlGaAs-pGaAs interface, the diffusion length in pGaAs, and the lifetime in 
pGaAs were critical parameters which help determine the spectral response 
at the blue end of the spectrum. The parameters in the nGaAs layer which 
determine the spectral response in the red end of the spectrum were found 
to be the diffusion length, diffusion coefficient, and lifetime. 

Annealing models were constructed assuming both monomolecular and 
bimolecular reactions for the interstitial-vacancy recombination reactions. 
These theoretical models produced curves which agree quite well with the . 
experimental data. 
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